Snakemake Report

Loading Snakemake Report...

Please enable Javascript in your browser to see this report.

Loading 3.2 MB. For large reports, this can take a while.

Workflow

human_genomics_pipeline processes paired-end sequencing data (WGS or WES) using bwa, sambamba and GATK4. It is based on the GATK4 best practices for data pre-processing for variant discovery

Run parameters:

Input data type: Single
Reference genome: /store/lkemp/publicData/referenceGenome/gatkBundle/GRCh37/ucsc.hg19.fasta
dbSNP database /store/lkemp/publicData/dbSNP/gatkBundle/GRCh37/dbsnp_138.hg19.vcf

Click the nodes to obtain details about each step.

Base recalibration

Quality checks

Trimming

Statistics

If the workflow has been executed in cluster/cloud, runtimes include the waiting time in the queue.

Jobs 1

Input files

qc/fastqc/NIST7035_NIST_chr1_R1_fastqc.zip
qc/fastqc/NIST7035_NIST_chr1_R2_fastqc.zip

Output files

qc/multiqc/pre_trim_multiqc_report.html

Conda software stack

bioconda::multiqc =1.8

Code

1	`( time multiqc {input} -o qc/multiqc/ -i pre_trim 2> multiqc_pre_trim.stderr ) 2> times/multiqc_pre_trim_time.txt`

Jobs 1

Input files

trim_galore/NIST7035_NIST_chr1_R1_val_1_fastqc.zip
trim_galore/NIST7035_NIST_chr1_R2_val_2_fastqc.zip

Output files

qc/multiqc/post_trim_multiqc_report.html

Conda software stack

bioconda::multiqc =1.8

Code

1	`( time multiqc {input} -o qc/multiqc/ -i post_trim 2> multiqc_post_trim.stderr ) 2> times/multiqc_post_trim_time.txt`

Jobs 1

Input files

mapped/{sample}_bwa_recal.bam
/store/lkemp/publicData/referenceGenome/gatkBundle/GRCh37/ucsc.hg19.fasta
/store/lkemp/publicData/dbSNP/gatkBundle/GRCh37/dbsnp_138.hg19.vcf

Output files

vcf/{sample}_raw_snps_indels_AS_g.vcf

Conda software stack

bioconda::gatk4 =4.1.6.0

Code

( time gatk HaplotypeCaller -I {input.bams} -R {input.refgenome} -D {input.dbsnp} -O {output} --tmp-dir {params.tdir} --native-pair-hmm-threads {threads} {params.padding} {params.intervals} 2> gatk4_HaplotypeCaller_single.stderr ) 2> times/gatk4_HaplotypeCaller_single_time.txt

Jobs 1

Input files

../fastq/{sample}_R1.fastq.gz
../fastq/{sample}_R2.fastq.gz

Output files

qc/fastqc/{sample}_R1_fastqc.html
qc/fastqc/{sample}_R2_fastqc.html
qc/fastqc/{sample}_R1_fastqc.zip
qc/fastqc/{sample}_R2_fastqc.zip

Conda software stack

bioconda::fastqc =0.11.9

Code

1	`( time fastqc {input.R1} {input.R2} --outdir qc/fastqc 2> {log} --threads {threads} 2> fastqc.stderr ) 2> times/fastqc_time.txt`

Jobs 1

Input files

../fastq/{sample}_R1.fastq.gz
../fastq/{sample}_R2.fastq.gz

Output files

trim_galore/{sample}_R1_val_1.fq.gz
trim_galore/{sample}_R2_val_2.fq.gz
trim_galore/{sample}_R1_val_1_fastqc.zip
trim_galore/{sample}_R2_val_2_fastqc.zip
trim_galore/{sample}_R1.fastq.gz_trimming_report.txt
trim_galore/{sample}_R2.fastq.gz_trimming_report.txt

Conda software stack

bioconda::trim-galore =0.6.5

Code

1	`( time trim_galore --paired {input.R1} {input.R2} -o trim_galore/ --cores {threads} {params} 2> trim_galore_pe.stderr ) 2> times/trim_galore_pe_time.txt`

Jobs 1

Input files

mapped/{sample}_sorted_mkdups_rgreplaced.bam
mapped/{sample}_sorted_mkdups_rgreplaced.bam.bai
mapped/{sample}_recalibration_report.grp
/store/lkemp/publicData/referenceGenome/gatkBundle/GRCh37/ucsc.hg19.fasta

Output files

mapped/{sample}_bwa_recal.bam

Conda software stack

bioconda::gatk4 =4.1.6.0

Code

1	`( time gatk ApplyBQSR -I {input.bams} -bqsr {input.recal} -R {input.refgenome} -O {output} {params.padding} {params.intervals} 2> gatk4_ApplyBQSR.stderr ) 2> times/gatk4_ApplyBQSR_time.txt`

Jobs 1

Input files

mapped/{sample}_bwamem_sorted_mkdups.bam

Output files

mapped/{sample}_sorted_mkdups_rgreplaced.bam

Conda software stack

bioconda::gatk4 =4.1.6.0

Code

1	`( time gatk AddOrReplaceReadGroups -I {input} -O {output} --TMP_DIR {params.tdir} {params.other} 2> gatk4_AddOrReplaceReadGroups.stderr ) 2> times/gatk4_AddOrReplaceReadGroups_time.txt`

Jobs 1

Input files

mapped/{sample}_sorted_mkdups_rgreplaced.bam

Output files

mapped/{sample}_sorted_mkdups_rgreplaced.bam.bai

Conda software stack

bioconda::sambamba =0.7.0

Code

1	`( time sambamba index -p {input} -t {threads} 2> sambamba_index_rgadd.stderr ) 2> times/sambamba_index_rgadd_time.txt`

Jobs 1

Input files

mapped/{sample}_sorted_mkdups_rgreplaced.bam
mapped/{sample}_sorted_mkdups_rgreplaced.bam.bai
/store/lkemp/publicData/referenceGenome/gatkBundle/GRCh37/ucsc.hg19.fasta
/store/lkemp/publicData/dbSNP/gatkBundle/GRCh37/dbsnp_138.hg19.vcf

Output files

mapped/{sample}_recalibration_report.grp

Conda software stack

bioconda::gatk4 =4.1.6.0

Code

( time gatk BaseRecalibrator -I {input.bams} -R {input.refgenome} --known-sites {input.dbsnp} -O {output} {params.padding} {params.intervals} 2> gatk4_BaseRecalibrator.stderr ) 2> times/gatk4_BaseRecalibrator_time.txt

Jobs 1

Input files

mapped/{sample}_bwamem_sorted.bam

Output files

mapped/{sample}_bwamem_sorted_mkdups.bam
mapped/{sample}_bwamem_sorted_mkdups.bam.bai

Conda software stack

bioconda::sambamba =0.7.0

Code

1	`( time sambamba markdup -p {input} {output.bam} --tmpdir={params.tdir} {params.other} -t {threads} 2> sambamba_mkdups.stderr ) 2> times/sambamba_mkdups_time.txt`

Jobs 1

Input files

mapped/{sample}_bwamem.bam

Output files

mapped/{sample}_bwamem_sorted.bam
mapped/{sample}_bwamem_sorted.bam.bai

Conda software stack

bioconda::sambamba =0.7.0

Code

1	`( time sambamba sort -p {input} -o {output.bams} --tmpdir={params.tdir} {params.other} -t {threads} 2> sambamba_sort.stderr ) 2> times/sambamba_sort_time.txt`

Jobs 1

Input files

/store/lkemp/publicData/referenceGenome/gatkBundle/GRCh37/ucsc.hg19.fasta
trim_galore/{sample}_R1_val_1.fq.gz
trim_galore/{sample}_R2_val_2.fq.gz

Output files

mapped/{sample}_bwamem.bam

Conda software stack

bioconda::samtools =1.10
bioconda::bwa =0.7.17

Code

1	`( time bwa mem {input.refgenome} {input.R1} {input.R2} -t {threads} -M \| samtools view -@ {threads} -Sbh - > {output} 2> bwa_map.stderr ) 2> times/bwa_map_time.txt`